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The QCD treatment of the photon structure beyond the leading order is discussed with emphasis on a proper choice of the 
factorization scheme and/or the boundary conditions. Recent parametrizations of the photon's parton content are examined 
and compared. The sensitivity of the photon structure function on the QCD scale parameter is reconsidered. 
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1. Introduction 

The hadronic structure of real photons [|l|, see 
also 0, is most easily illustrated by considering 
inclusive electron-photon deep-inelastic scattering 
(DIS). This process is completely analogous to the 
usual lepton-nucleon DIS. Consequently, the unpo- 
larized cross section for 67 — > e"i*"f -^ eX can 
be expressed in terms of two structure functions 
F^'^iXjQ'^), where x denotes the Bjorken scaling 
variable and Q^ is the virtuality of the probing off- 
shell photon. Only F2 has been measured so far. 
The main issues addressed in this talk can already be 
introduced at the parton model level. Here, two dif- 
ferent contributions occur, which originate in point- 
like and non-pointlike photon-quark couplings. 

The pointlike contribution derives from the lowest 
order QED process ('Born-Box') 7*7 ^ (79 0, and 
is for Q'^ 3> m^ given by 



-^2T.-£E24^ 



(0) 



X 



ia:)ln^ + C,{x)] . (1) 



Effective masses ruq of the quarks with charges e^ 
had to be introduced in order to avoid coUinear sin- 
gularities. Contrary to all hadrons the x-dependence 
of F2, governed by kg '{x) and C-y{x), is calculable 
and the structure function rises linearly with In Q^ . 
Analogous to the usual MS treatment of hadronic 
DIS, the quark densities in the photon can be de- 
fined by absorbing the universal mass singularity, 

^Il« = ^llTs = £^^3[.^-K(l-.nin^,(2) 

where the concrete form of kg has been put in. The 
Cy-term, not absorbed into the parton distributions 



in (2), then acts as the subleading, 'direct' contribu- 
tion to F2 ■ However, this separation creates prob- 
lems in the next-to-leading order QCD treatment, 
mainly since C-y is negative and divergent at a; — > 1. 
Solutions to this problem, either by introducing a 
different factorization scheme ('DIS^') or by adding 
some 'technical' MS input for the evolution equa- 
tions, are examined in section 2. 

The hadronic, non-pointlike contribution to F^ 
can be illustrated by the well-known coupling of the 
photon to the vector mesons p, to and (j) (vector me- 
son dominance, VMD) [n], resulting in 



^2,VMD 



= (4^a//2)K 



(3) 



This part is completely analogous to the behaviour 
of hadrons: there is no rise with In Q^ and the x- 
shape is not calculable in perturbation theory. Hence 
unknown input parton distributions enter the de- 
scription just as in hadron structure. In section 3 re- 
cent parametrizations of q'^ and the less well known 
gluon density G^ are discussed and compared. 

The dependence of F2 on the unknown effec- 
tive quark masses is removed by the inclusion of 
the dominant QCD corrections, leading to the fa- 
mous large-x large-Q^ parameter-free asymptotic 
QCD predictions [00]. These results, however, are 
not appropriate at scales accessible at present or in 
the near future, Q^ < 300 GeV^ g. The original 
hope on an especially clean as determination from 
F2 has therefore been dampend. Some debate has 
arisen on how much sensitivity survives in a proper, 
non-asymptotic treatment including the uncalcula- 
ble hadronic boundary conditions, see e.g. 0. A 
short note is added to this discussion in section 4. 
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2. QCD: Next to Leading Order 



The generalized evolution equations for the 0{a) 
quark and giuon content of the photon are |g|,^ 



dlnQ2 



27r 
27r 



as 

27r 



[^< 



GG ' 



G' 



(4) 



Here denotes the usual convolution, and the nor- 
mal hadronic part has not been written out com- 
pletely. The photon distribution in the photon, F^ = 
(5(1 — x) — 0{a), has dropped out at the order of the 
electromagnetic coupling a considered here, result- 
ing in the characteristic inhomogeneous fc-terms. A 
salient consequence of (4) is the absence of a mo- 
mentum sum rule interrelating q"^ and G '^ . Hence an 
important constraint on the gluon density, present 
in the hadronic case, is missing here. The splitting 
functions Pij and ki are up to now known to next-to- 
leading order (NLO) accuracy in the strong coupling 
as = as{Q^), see [|| and ]lO| , p^ , respectively. 
The singlet distributions can be decomposed as 



q'' = 






— <lpL + ^had ' 



(5) 



where q"' =q'^ has been used. The well-known homo- 
geneous hadronic solution Q/j^^ contains the pertur- 
batively uncalculable boundary conditions q '''(Qq). 
The inhomogeneous ('pointlike', PL) solution then 
satisfies qpi^iQo) = 0. 

An analytical solution of (4) is possible for Mellin- 
moments. The NLO pointlike part reads |ll| , p2[ 

-{- + u}{^-[as/as{Ql)r'}-^y 



QPL 



: a 



{ 



l-ias/asiQDYl-.b + Oias) , (6) 
J a 



with a, 6, d and U being known Mellin-?i depen- 
dent combinations of the splitting functions. The 
cc-dependent distributions are obtained by a numer- 
ical Mellin-inversion of (6) and the corresponding 
hadronic and non-singlet expressions. A major ad- 
vantage of this procedure is that inconsistent higher 
order 0{as) terms, which would lead to spurious 
a; — > 1 singularities, can be trivially omitted. In a 
direct NLO x-space calculation, their cancellation 
requires somewhat cumbersome algorithms pSlMl . 



- 2, PL 

NLO for the MS and DlS-y factorization schemes. Also shown 

is the modified F^ pj^, of jisj (AFG). Qq = 1 GeV^, three 

active flavours and A = 0.2 GeV have been used. 

The structure function F2 is at NLO given by 



V = V2e^L^-H^(G,®g^ 



Cg(E>G^) 






(7) 



where the summation extends over the light u, d 
and s quarks. The heavy flavour contribution F^ f^ 
has recently been calculated to second order in as 
[ p^ , see also [|l^. Cq^a are the usual hadronic co- 
efficient functions, and for the commonly used MS 



factorization one has G-, = 6 G^: 



The latter 'di- 



rect' photon term, however, causes difficulties in this 
scheme, since it leads to a large LO/NLO difference 
for the pointlike contribution F^ pj^ shown in Fig. 
1. Especially, it is negative and divergent at large x: 

C^j^^?,[\n{l-x)-l] for x^l . (8) 

This apparent perturbative instability in F2 has to 
be compensated by the NLO MS input parton den- 
sities, which therefore are forced to be very differ- 
ent from the LO ones. In particular, a physically 
motivated VMD input shape, vanishing for a; ^ 1, 
cannot be employed in NLO here. 
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These problems are avoided by working in the 
DIS^ factorization scheme introduced in fl^. Here, 
the Cj term is (re-)absorbed in the quark densities, 
i.e. the constant and the logarithmic parts of (1) are 
not, as in (2), separated for the definition oi q'' . The 
connection to the usual MS scheme is 



7 

Idis-, 






'MS 



C 



-l,DIS^, 







(9) 



and the coefficient functions Cq and Cg as well 
as the definition of the gluon density remain un- 
changed, in contrast to the hadronic DIS scheme 
[ pj| . Thus in DIS^ F^ retains the usual hadronic MS 
form without a direct term also at NLO, resulting in 
a good LO/NLO stabihty of F^p^ (Fig. 1). Hence 
in this scheme perturbatively stable, physically mo- 
tivated inputs for the photonic parton densities, e.g. 
from VMD as in ^^^, can be used in NLO. 

Due to the non-universality of the coefficient func- 
tion C-y, F2 has been assigned a special role in the 
redefinition (9) of the quark distributions. The first 
(pragmatic) motivation for doing so is that F2 is the 
only photonic DIS structure function measured so 
far and in the near future. The second (theoretical) 
reason stems from the analogy with the hadronic 
DIS scheme, where F2 is special because of its con- 
nection with the charge sum rules |17[ . 

An additional advantage of the DIS^ scheme, im- 
proving the perturbative stability of the evolution, 
is that the strongest MS poles for x -^ in the NLO 
photon-parton splitting functions k\ disappear: 

DIS^, 
(10) 
A corresponding cancellation takes place in the anal- 
ogous timelike case, where it is even more important: 
besides the corresponding structure functions, also 
the photonic gluon fragmentation function D^ pj^ is 
strongly negative at small and medium x in the MS 



k^ 





scheme |20 . These problems are removed by using 
the timelike DIS7 factorization [ pl[ . 

An equivalent MS formulation of the above solu- 
tion to the C-y-problem is of course possible and has 
been approximately employed in p2[ . It can be writ- 
ten as a modification ('PL') of the pointlike part in 
(5) by an additional 'technical' NLO quark input. 



qy^iQi) 






(11) 



NLO for the DIS^ scheme compared to the physically equiv- 
alent MS distribution (PL) with the input (11). The param- 
eters Qq etc. are as in Fig. 1. 

This leads to F'^ — (Qo) = ^ and thus allows for 
similar 'physical' inputs to q^^^ at LO and NLO. 
The resulting quark densities, however, depicted in 
Fig. 2 together with their pointlike LO and DIS^ 
counterparts, exhibit a rather unphysical shape. 

In a consistent NLO calculation, where e.g. the 
input (11) is evolved and convoluted only with LO 
quantities, the above MS procedure is strictly equiv- 
alent to the DIS^ treatment for all physical quanti- 
ties. On the other hand, as soon as not all higher 
order 0{as) terms are discarded, the MS factoriza- 
tion turns out to be much less well-behaved than 
the DIS^ scheme. This is obvious from Fig. 3, where 
the results of a 'slightly' inconsistent calculation, in 
which only the additional 0{as) terms arising from 
the convolutions like asCq^q'' in (7) have not been 
omitted, are compared to the consistent result for 
both schemes. Consequently, in particular for appli- 
cations in the MS scheme, parametrizations of NLO 
photonic parton densities should, as in [|8|, supply 
also the leading part of q^^^o ^"^ (^) (v^ 9lo) neces- 
sary for consistent calculations. 



A. Vogt / Photon structure: QCD treatment and parton densities 



3. Photonic Parton Densities 



Fig. 3. The consistent FS' pr of Fig. 1 in comparison with 
inconsistent calculations, where spurious 0{as) terms from 
the convolutions in (7) have not been cancelled. 

Very recently, an alternative approach has been 
suggested |g^, where a different (universal) techni- 
cal input has been constructed by a detailed analysis 
of the momentum integration of the Box-diagram: 



QlL'iQl) 






Gl^,iQl)^0 



(12) 



with Ci^ 



3 ln(l — x) for X -^ 1. The resulting mod- 
ified pointlike structure function F^ pj^, , also pre- 
sented in Fig. 1, remains negative at large x due to 
the uncompensated -1 in (8) only very close to Qq. 
However, a significant dissimilarity to the LO result 
remains, which demands compensation by sizeably 
different LO/NLO hadronic inputs. 

Finally, it should be mentioned that the theoret- 
ical uncertainty of NLO photon structure calcula- 
tions has been investigated recently [p5|J23J. It has 
been found that a variation of the renormalization 
and factorization scale M^ from Q^/A to AQ^ leads 
to an increase of F^ by less than 15% (80%) and 
5% (40%) at Q2 ^ 6 GeV^ and Q^ = 50 GeV^ in 
NLO (LO), respectively [Q. Thus theory is ready 
for significant comparisons with future high-Q^ pre- 
cise data, e.g. from LEP 200. 



In order to specify the photon's parton distri- 
butions, one has to choose a reference scale Qq 
and to fix the perturbatively uncalculable bound- 
ary conditions q^iQ^), G'^{Qq) for the evolution 
equations (4). Experimentally, q"' (mainly u^) is 
known at 0.05 ~ a: ~ 0.8 from F^ measurements 
[||,|5|J2^,|2^,||,|9|||l with an accuracy of about 
10 -^ 20 %. On the other hand, little is known about 
C up to now: there is evidence for C^ ^ 0, and 
an extremely huge and hard gluon (LAC 3, see be- 
low) has been ruled out by jet production data 
pl| , ^2|j3^ , |3^ , see also ||3|] . The prospects on a mea- 
surement of the gluon density are discussed in [Q . 
This lack of experimental information is especially 
serious here since, as explained above, there is no 
energy- momentum sum rule relating q'^ and C . 

VMD provides a connection between the quark 
and gluon inputs by the hadronic momentum sum 
rule. Moreover, by SU(6) arguments, G''{Q'^) can 
then be estimated from the pionic gluon density 
C" , which has been determined from direct pho- 
ton production 7rp -^ ^X ||37| , |38| . However, for the 
traditionally chosen input scales, Q^ > 1 GeV^, a 
pure VMD ansatz is known to be insufficient. An 
additional hard quark component has to be supple- 
mented in order to achieve agreement with the F2 
data at larger Q^ ^Q- In contrast to the quarks, 
where the Born-Box (1) can provide a natural addi- 
tional input, a corresponding natural choice beyond 
VDM does not exist for the gluon density. 

In view of this situation two approaches have 
been used. First, one can keep Qo > \ GeV, fit 
the quark densities to F2 data and 'guess' the 
gluon input. This method has been adopted in |39[ 
and, more recently, in ^,Q. The second possi- 
bility is to try to maintain the VMD idea and to 
start the perturbative evolution at a very low scale 
Qo < 1 GeV [^J23[|^,||. In the following, the re- 
cent parametrizations are discussed. All of them use 
a QCD scale parameter of A = 200 MeV for four 
active flavours in LO and NLO. The PEP, PETRA 
and TRISTAN data on F2 available when most of 
these parametrizations were done are shown in Fig. 
4. Since then, also first results from LEP have been 
pubhshed M, see also l43|. The quark and gluon 
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Fig. 4. The 1991 world data on F^ p3,E|,E|J2^J28lp|] compared to the results obtained from the LO and NLO GRV parametriza- 
tions, where the charm contribution has been calculated via the lowest order Bethe-Heitler process llSl. The controversial [M 
low-Q^ large-x TPC/27 data points close to the resonance region (open circles) have not been used for the determination of k in 
(13). 



densities are compared in Fig. 5 and 6, respectively. 

LAC (only LO) @: 

The input scales are Q§ = 4 GeV for sets 1, 2 and 
Qq = 1 GeV^ for set 3. The boundary conditions 
have been fitted to the data of Fig. 4 at Q^ > Qq 
without any recourse to physical constraints on the 
quark flavour decomposition and on the gluon den- 
sity. This has partially led to rather unphysical re- 
sults even in the quark sector, e.g. u'^ — l.2a, 



dT = 0.95a, s^ 



2.8a and 



1.6a for LACl at 



a: = 0.1 and Qq — 10 GeV . Moreover, such a pro- 
cedure gives rise to wild reactions of the fitted gluon 
density on fluctuations and offsets in the data: in 
LACl and LAC 2, Gt is very large at x < 0.1 due 



to the small- a; interplay of PLUTO data 1 24 at low 
Q2 with JADE H and prehminary CELLO || re- 
sults at larger Q^ in the fit. Likewise, the extremely 
huge and hard LAC 3 gluon is driven by the inclusion 
of the large-a; small-Q^ TPC/27 data [|7|, which fie 
close to the resonance region and are therefore sub- 
ject to potentially large experimental correction fac- 



tors and higher twist contributions |Q. In any case 
the LAC gluons clearly illustrate that present F2 
data do not provide useful constraints on C*". 



GRV (LO/NLO) (18): 

The gluon and sea quark densities of the proton 
and the pion have been generated from a valence-like 
structure at a common, very low resolution scale, 
^i ~ 0.55 (0.5) GeV in NLO (LO) ^^. First, 



this procedure has the advantage that few param- 
eters in the gluon/sea sector are needed. The pionic 
gluon density, being a parameter-free prediction in 
the simplest version of this approach, shows very 
good agreement with the results from direct photon 
production [ p7|j3^ . Secondly, predictions at small x 
arise, where the distributions are determined by the 
QCD dynamics and do virtually not depend on am- 
biguous input assumptions. HERA results on F2 
M,E8| exhibit excellent agreement with these pre- 
dictions, especially when the charm contribution is 
properly taken into account |^ . 

Therefore it is natural to utilize these low scales 
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Fig. 5. Photonic M-quark densities at LO and NLO. The LO distributions |18lba,kcj] have been fitted to 1991 F^ data, for the 
different assumptions and data sets used see the discussions in the text. The NLO results [ |l8|J22y23| are presented in the MS 
scheme. Here only the GRV distribution has been directly fitted to data. 



Fig. 6. Photonic gluon distributions in LO and NLO. The spread of the LO curves ||lq,p2,M exaggerates the present experimental 
uncertainty on C", since LAC 3 has been ruled out [p]lj32|. On the other hand, the similarity of the NLO results [ |l8y22[|23| is due 
to common VMD prejudices and not enforced by data. 
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with as/TT w 1/4 also in the photon case for a pure 
VMD ansatz in LO and NLO (DIS^): 



iq\G''){^l') 



(13) 



K{ATra/fp){qj,o,G^o){fi ) , 

with the couphng /p/(47r) = 2.2 and 1 ~ k ~ 2 §. 
In (13) the pion distributions of |Q have been used 
instead of the experimentally unknown p densities. 
K accounts in a simple way for the higher mass vec- 
tor mesons p, to and </>. Since this parameter is not 
sufhciently specified by VMD, it has been fitted to 
the F2 data in Fig. 4. Good fits are obtained for 
the very reasonable values k = 1.6 (2.0) in NLO 
(LO). Data does not favour neither an additional 
hard input at /i^ nor a later onset of the pointlike 
contribution. Thus a pure VMD input at Qq = p? is 
in fact sucessful and allows for an approximate pre- 
diction of C" in this 'dynamical' approach. Uncer- 
tainties stem from the p — > tt substitution and from 
the pion-valence normalization [p8| , p0| entering the 
determination of k. While the former is difficult to 
quantify, the latter amounts to about 20%. 

It is interesting to note that the same scale Qo = 
0.5 GeV for a LO VMD input has been obtained 
from quite different 7p total cross section considera- 
tions |42|. Consequently, the resulting distributions 
resemble the ones of [Q . In order to remain flexible 
on the VMD part, the LO pointlike solution (5) has 



been parametrized seperately in |42 



GS (LO/NLO) ||: 

The input scale is Qq = 5.3 GeV^, the average of 
the 1ow-Q2 PLUTO data set ^. The F^'' measure- 
ments at Q^ > Qq in Fig. 4 have partially been used 
to fit the free parameters indicated below. In view 
of the gluon ambiguity discussed above, two sets are 
provided with different assumptions on G'*'. At LO, 
the inputs of sets 1 and 2 are parametrized as 



q^2 = VMD + QPM(m,) 



(14) 



G^ = VMD + —Pg^°^ (g)J:^ ,G^ ^ VMD 



The VMD part is treated as in (13), with the old 
pion distributions from Q being used. The effective 
quark masses ruu^d.s in the parton model (QPM) ex- 
pression (1) and K in (13) are determined from data. 
The difference between G^ and G2 presented in Fig. 
5 is about twice as big as the LO pointlike contri- 



bution evolved from the start scale Qq — 0.5 GeV 
employed in ||l|,|2). Corresponding NLO (MS) den- 
sities have been constructed without a new fit by 
enforcing the same i^2^(Qo) ^^ i^- LO together with 
assumptions on the flavour decomposition. 

AFG (only NLO) |2|: 

These distributions supersede the NLO ones of 
|ll| , where a pure VMD input was employed at LO 
and NLO (MS) at Qo = 0.5 GeV. Here, the techni- 
cal input (12) is supplemented by a VMD ansatz at 
Qo ~ 0.7 GeV « nip. No fit to data is performed. 
Somewhat different from (13), a coherent sum of p, 
uj and 0, is used, where identical valence and gluon 
distribution in all three vector mesons and the pion 
have been assumed: 



u'* = Ka 



9 



^val 



3 



G'^=Ka-G'' (15) 



The pion densities are adopted from [Q with the 
normalization K left free in the parametrization. 
The standard choice reads K ^ I, with a 50 % vari- 
ation allowed by high-Q^ data on F2 ■ 



4. as and the photon structure 

The sensitivity of F2 to the QCD scale param- 
eter A considerably decreases with increasing scale 
Qq, at which the evolution (4) of the photonic par- 
tons is starlet H,B3,p3[. For Q^ ~ 5 GeV^, as in 
[P4I , only a marginal effect is left. It interesting, how- 
ever, to reconsider this issue employing the 'dynami- 
cal' approach, where the parton densities of nucleon, 
pion and photon are related by valence-like inputs 
and VMD at a common, very low resolution scale, 
p ~ 0.55 GeV for A^IL = 0.2 GeV |l|,||,|6|. There- 
fore, this procedure, carried out only for this fixed 
A so far, has been repeated for different values of 
Ap-g with /i(A) fixed by proton structure informa- 
tion. The free parameter in the photon sector, k in 
(13), has been fitted to the pubhshed results on F2 
with a cut of Q^ > 4 GeV^ in order to suppress pos- 
sible higher- twist contributions. These fits exhibit a 
rather good A-sensitivity in this dynamical scenario, 



A(l^=(210±50e.p.)MeV 



(16) 
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where the 1 a error given in (16) arises from the total 
experimental uncertainties. A more detailed analysis 
including theoretical uncertainties will be presented 
in a seperate paper ||54| . Since more precise data are 
expected, especially from LEP 200, prospects do not 
seem to be too dim for a competitive as determina- 
tion from F2 ■ 
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